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Abstract. We consider the determination of a s from r hadronic decays, by investigating the contour- 
improved (CI) and the fixed-order (FO) renormalization group summations in the frame of a new per- 
turbation expansion of QCD, which incorporates in a systematic way the available information about the 
divergent character of the series. The new expansion functions, which replace the powers of the coupling, 
are defined by the analytic continuation in the Borel complex plane, achieved through an optimal con- 
formal mapping. Using a physical model recently discussed by Beneke and Jamin, we show that the new 
CIPT approaches the true results with great precision when the perturbative order is increased, while the 
new FOPT gives a less accurate description in the regions where the imaginary logarithms present in the 
expansion of the running coupling are large. With the new expansions, the discrepancy of 0.024 in a s (m 2 ) 
between the standard CI and FO summations is reduced to only 0.009. From the new CIPT we predict 
(m 2 ) = 0.320 

-o'oo9> which practically coincides with the result of the standard FOPT, but has a more 
solid theoretical basis. 



1 Introduction 

The precise determination of a s from r hadronic decays is 
one of the most important results in perturbative QCD pQ . 
The subject has been treated by many authors (see [2]-[T2] 
and references therein). Among the most important im- 
provements we mention the so-called "contour-improved 
renormalization-group summation" [4], [5], which avoids 
the large logarithms of the usual "fixed-order" expansion 
along the integration contour relevant for a s extraction. 
The higher orders of perturbation theory, expressed by the 
so-called renormalons, were also investigated, especially as 
concerns their effect on the precision of the theoretical de- 
terminations [7j, 0. 

The problem was revisited recently, after the calcula- 
tion of the Adler function up to fourth order [T3] , the same 
at which the function describing the running of the cou- 
pling is known [14j , [15] . The determination of a s from the 
ALEPH spectral function data was reconsidered in [IT)] , 
where the prediction based on CIPT a s (m 2 ) = 0.344 ± 
0.009 was made. On the other hand, an updated version of 
FOPT [T7] led to the prediction a s {m 2 r ) = 0.320±£oo7- As 
remarked in [T7], the discrepancy of 0.024 between CIPT 
and FOPT appears to be the largest systematic theoreti- 
cal uncertainty in the a s determination, and it does not go 
away by adding the presently known higher-order terms. 

The two summation methods were analyzed in de- 
tail in [T7] by means of a "physical" model for the Borel 



transform of the Adler function. The conclusion of this 
study was that, somewhat surprisingly, FOPT is prefer- 
able, CIPT failing to approach the true result, although 
it is a priori more consistent. 

The purpose of the present work is to understand bet- 
ter the relation between the CI and FO summation meth- 
ods and their consequences on the extraction of a s . To this 
end we investigate these two methods using a new per- 
turbation expansion advocated by us some time ago [18] . 
The new expansion separates the intrinsic ambiguity of 
the perturbation theory due to the infrared regions of the 
Feynman diagrams, from the divergent character of the se- 
ries. This is achieved if the usual perturbative expansion 
is replaced by a series with better convergence properties 
(each new term improving the accuracy of the approxi- 
mation instead of spoiling it). As shown in [IB]- [ID], the 
new expansion is defined by using the analytic continua- 
tion in the Borel complex plane by an optimal conformal 
mapping. The properties of the novel expansion functions 
were analyzed in detail in |19j.|20j and some applications 
were discussed in [18], [2T] , 

The paper is organized as follows: in Section [2] we 
briefly review the CI and FO summations of the Adler 
function relevant for r decays. In Section [3] we discuss the 
Borel transform and in Section [4] we present the CI and 
FO versions of the new perturbation theory for the Adler 
function. In Section [5] we apply these two methods to cal- 
culate a s (m 2 ), and show that the difference between the 
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two predictions is significantly smaller than in the stan- 
dard case. In order to understand these results, in Section 
[6] we investigate the new expansions using as an example 
the physical model proposed in [17]. More comments on 
the method of conformal mapping and its relevance for 
the physical case are made in Section [7] We summarize 
our results in Section [8] where we emphasize that the CI 
summation combined with the new perturbation expan- 
sion leads to a precise determination of a s (m^). 



2 Standard CIPT and FOPT 

The relevant quantity for the extraction of a s (m^) is the 
integral 

*(°) = -L I *L u (s)D(b), (1) 
2m J s 

\s\=s 

where sq = m^, w(s) = 1 — 2s/s + 2(s/sq) 3 — (s/sq) 4 and 
D(s) — D(s) — 1 is the reduced Adler function in massless 
QCD. It is written formally as the renormalization-group 
improved serie^ 



D(s) 



OO 

E 

n=l 



K n (a s (s)Y 



(2) 



where a s (s) = a s (s)/7r. In the MS scheme, for nf = 3, 
the coefficients K n calculated up to now have the values: 

Ki = 1, K 2 = 1.6398, K 3 = 6.3712, K± = 49.076. (3) 

Several methods were proposed for estimating the higher- 
order perturbative coefficients from the low-order ones 
[22], [23]. In particular, in [17] the authors adopt the value 



Kb = 283. 



(4) 



We mention that at large n the coefficients display a facto- 
rial increase, K n ~ n\, so that the series in |J5J) is divergent. 
In writing ^ we follow the convention often adopted in 
physical papers, writing the sign of equality even if the se- 
ries on the right hand side is divergent and the equality is 
impossible. According to Dyson's proposal [53] from 1952, 
the series is then regarded as asymptotic to D(s) standing 
on the left hand side. If the series is convergent, the rela- 
tion ([2]) is understood as equality. Analogous series in the 
text below (see ([16]), (|20]l, ([27]). ([30]) . etc.) are understood 
in a similar sense. 

The contour-improved (CI) summation amounts to in- 
troducing ([2]) in (fTJ) , and performing the integral with 
a s (s) calculated locally from the solution of the renormal- 
ization group equation. The solution is known at present 
up to four loops [H] , the first coefficients /3j of the (3 func- 
tion, calculated in the MS scheme and n/ = 3 being 

/3 = 9/4, 0X = 4, 2 = 10.0599, /3 3 = 47.228. (5) 



1 The normalization of D is that adopted in [T7], where K n 
are denoted as c„,i. For simplicity, in (|2]| the scale was set to 
£ = 1. The general case will be discussed at the end of Section 

H 



The Taylor expansion of a s {s) in terms of a reference 
point si reads 



J">1 



(6) 



where the coefficients £j depend on 771 = ln(s/si) and 
on the coefficients f3j. By inserting this expansion with 
si = sq in ([2]) and rearranging the expansion in powers 
of a s (so), one obtains the FO summation, which can be 
written as 



D(s) = J2[K n + K n (s)\ (a s (s )) n , 



(7) 



where K n (s) are polynomials of 770 = m (s/so), and depend 
on the coefficients 0j and Kj for j < n. 



3 Borel transform 

The series ([2]) can be formally written as the Borcl-Laplace 
transform 



D(s) = \- [ e -»/0ao».(-)) B(u) du , (8) 
Po J 


where B{u) is the Borel transform of the Adler function, 
defined by the power series 



B{u) 



n=0 



b n u r ' 



(9) 



with b n related to the original perturbative coefficients 
appearing in @ by 



K 



n+l 



n > 0. 



(10) 



According to present knowledge, the function B(u) has 
branch point singularities in the w-plane, along the neg- 
ative axis - the ultraviolet (UV) renormalons - and the 
positive axis - the infrared (IR) renormalons [25]. Specifi- 
cally, the branch cuts are situated along the rays u < — 1 
and u > 2. The nature of the first branch points was es- 
tablished in [26] and in [27] (see also [17]). Thus, near the 
first branch points, i.e. for u ~ — 1 and u ~ 2, respectively, 
B{u] behaves as 



B{u) 



n 



(i + uy 



-, B{u) 



(1 - w/2) 72 



(11) 



where the residues T\ and r2 are not known, but the ex- 
ponents 71 and 72 are known [5S] ; [27] . 

Due to the singularities along the positive axis, the in- 
tegral |5]| does not exist. The ambiguity in the choice of 
prescription is often used as a measure of the uncertainty 
of the calculations in perturbative QCD. It is convenient 



I.Caprini and J. Fischer: a s from r decays 



3 



to define the integral by the Principal Value (PV) pre- 
scription: 



1 



£( s )ee-PV 
Po 



-u/(0 o a B ( S )) 



B(u) du . 



(12) 



where 



PV 



DC 

J f(u)du : 



- / [f(u + ie)du + f(u-ie)du] 



As discussed in [28], [29], the PV prescription is the best 
choice if one wants to preserve as much as possible the an- 
alyticity properties of the correlators in the s-plane, which 
are connected with causality and unitarity. 



4 New CIPT and FOPT 

In order to define a new perturbative expansion of the 
Adler function we shall apply the method of conformal 
mapping [30] . This method is not applicable to the series 
([2]), because D(s) (regarded as a function of a s (s)) is sin- 
gular at the point of expansion a s (s) = 0. The method 
can, on the other hand, be applied to ©, because B(u) is 
holomorphic at u = 0. 

We note that the expansion ([9]) converges only in the 
disk |u| < 1. A series with a larger domain of convergence 
can be obtained by expanding B(u) in powers of a new 
variable. As shown in [30], the optimal variable coincides 
with the function that performs the conformal mapping 
of the whole analyticity domain of the expanded function 
onto a disk in the new complex plane0. 

To find the explicit form of the optimal conformal map- 
ping, one should know the location of all the singularities 
of the Borel transform B(u) in the complex Borel plane. 
Unfortunately, present evidence of these singularities is 
scarce: the known singularities (IR and UV renormalons 
and instanton-antiinstanton pairs) are produced only by 
a subclass of Feynman diagrams, while the effect of all 
the diagrams on the nature and location the singulari- 
ties is not known. In the lack of rigorous results, addi- 
tional assumptions are made or special models are built. 
As for additional assumptions, we point out that univer- 
sally adopted has been to assume that B{u) has only the 
above-mentioned singularities on the real axis with a gap, 
being holomorphic elsewhere. 

Under this assumption, the optimal variable defined in 
130] reads [18]: 



w(u) 



VTTu- yl -u/2 

VTTTl+ y/l -u/2 



(13) 



2 For QCD, the use of a conformal mapping in the Borel 
plane was suggested in i3F and was applied in a more limited 



context in 
also in 33 



[32]. 



Applications of the method were considered 



This function maps the u-plane cut for u > 2 and u < — 1 
onto the unit disk |w| < 1 in the complex plane w — w(u), 
such that w(0) — 0, w(2) — 1 and w(— 1) = — 1. It is useful 
to give also the inverse u = u(w) of ([13]) : 



u(w) 



8w 



(14) 



3 - 2w + 3w 2 
According to general arguments [30j . the expansion 

B{u) = Y. d ^ wn ( 15 ) 



n>0 



converges in the whole disk \w\ < 1. Moreover, as shown 
in [3D] , the expansion ([T5"]) has the best asymptotic rate of 
convergence compared to all the expansions of the function 
B(u) in powers of other variables. 

The series ([15]) can be used to define an alternative 
expansion of D(s). This is obtained formally by inserting 
(|15p into (|12p and interchanging the order of summation 
and integration. Thus, we adopt the modified CIPT ex- 
pansion defined as [l8]-[20] 



D(s) = J2d n W n (s), 



(16) 



n>0 



where 



WJs) = -!-PV (e- u/{f3oa ° {s)) w n du. 
Po J 





(17) 



We empahsize that the expansion (fT5|) exploits only the 
location of the singularities in the Borel plane. However, 
in our case some information exists also about the nature 
of the singularities. We note that Eq. (fTT]) expressed in 
the variable w = w(u) implies 



B(u) 



(8/3)^1 
(1 + w) 2 ^ ' 



B(u) 



(4/3)^r 2 
(1 - w) 2 ^ 



(18) 



for the behaviour near the points w = — 1 and w = 1, 
respectively. The expansion (fl~5]) is expected to describe 
these singularities if a large number of terms is used. How- 
ever, since the nature of the singularities is known, it is 
convenient to incorporate it explicitly. This is achieved, 
for instance, by expanding the product (1 + w) 2,7l (l — 
w) 212 B{u) in powers of the variable w: 



(1 + w) 271 (1 - w) 2 ~< 2 B(u) = C n 



(19) 



n>Q 



The expansion (|19p converges in the whole disk \w\ < 1, 
i.e. in the whole complex u-plane, up to the cuts along 
the real axis. Moreover, since the singular behaviour of 
B(u) at the first branch points is compensated by the first 
factors in ([T9| . the series is expected to converge faster 
than ([T5]) . Also, the behaviour near the first singularities 
holds even for truncated expansions, which are used in 
practice. 
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It is important to stress that, while the expansion (|15jl 
is unique, the explicit inclusion of the first singularities 
of B{u) contains some arbitrariness. The description of 
the singularities by multiplicative factors is a possibility, 
but is not a priori necessary. Moreover, the factors are not 
unique. For instance, the information on the nature of the 
singularities can be exploited by factors in the u varia 
An advantage of the choice (|T^|) is that the multiplicative 
factors remain finite at large u. We will make more com- 
ments on this in Section [Jj 

The expansion ([T9"]) suggests the definition of the new 
CIPT 

£>(s)=^c n W„(s) ) (20) 



n>0 



where the expansion functions are defined as 



w" 



A) A, J (1+ W )^(1_ W )2 72 




du , 
(21) 

with w = w(u) defined in (|13[) . 

We emphasize that the expansions (|16[) and ([2"0"]) repro- 
duce the coefficients K n of the usual expansion @ , when 
the functions (fTTj) and ([2"T]) are expanded in powers of 
the coupling. In fact, as shown in [20], the new expansion 
functions are formally represented by divergent series in 
powers of the coupling, much like the expanded correlator 
itself. 

To obtain the FO version of the new expansions, we 
start from ([JJ and define the corresponding Borcl trans- 
form 



B(u,s) = J2b n {s)u n , 



(22) 



n=0 



where 



K{S) = pf^l ' n ^°- (23) 

Then the Adler function admits the formal representation 

oo 

D(s) = i- [ e- u/( -P° a ° {sa)) B(u, s) du . (24) 
Po J 
o 

By comparing Eqs.([22]), ([23]) with ©, |10]), we can write 



B(u, s) = B(u) +B l (u,s), 



(25) 



where Bi(u,s) is generated by the second term in the 
coefficients (|2T>1) . It follows that the singularities of B(u) 
are present also in the function B(u,s), which may have 
in addition singularities from the second term in (|25j) . In 
what follows we shall exploit the singularities of B(u), 
which are known, by expanding B(u, s) in powers of the 
variable w defined above: 



B(u,s) 



n>0 



d n (s) w T 



(26) 



3 More precisely, in 05] the product (l+«) 71 (1-m/2)' Y2 B(m) 
was expanded in powers of it, while in [21] the same product 
was expanded in powers of w. 



This leads us to the definition of a modified FOPT, anal- 
ogous to the CIPT expansion (|T5|) : 



(27) 



n>0 



in terms of the functions 



1 



W n = — PV / e -u/(Poa s (s )) w n du _ ( 2 g) 

PO J 


As we discussed above, it is convenient to impose explicitly 
the behaviour (jlip . which is done by expanding 

(1 + wf 11 (1 - w) 2 ^ 2 B(u, s) = ^2 5 «( s ) w n . (29) 

n>0 

Then, using (j2~4"|) we define the new FO expansion: 



D(s) = Vc„(s)W n , 



n>0 



in terms of the functions 



1 



(30) 



du . 



n (3 J (1 + w )2 7l (1 _ w )2 72 



( 31 ) 

In the next Section we shall determine a s using the 
new CI and FO perturbation expansions ([2"0"1) , ([2"Tj) and 
(I30p . (|3ip . which include in an explicit way the nature of 
the first singularities of the Borel transform. More com- 
ments on them and the general expansions (fT6|) . (fT7|) and 
(H7J, ([2S]), which replace the standard CIPT and FOPT, 
respectively, will be made in Section [7j 

By inserting ([2"0"]) in the integral ([T]), with a s (s) cal- 
culated from RGE applied locally, we obtain a new CI 
perturbation expansion for S^°\ Likewise, by using in the 
integral |T]) the expansion ([50]) . we obtain the new FO 
perturbation expansion of S^°\ 

We end this Section with a comment about the renor- 
malization scale. The starting point of our derivation was 
the renormalization group improved expansion which 
corresponds to the choice of the scale £ = 1 (with the 
notation used in [H]). The general case 



oo 



D(s) = yK n (0 (a s my 



(32) 



with the coefficients K n (0 given in [TT], is easily obtained 
by reordering ([2]) as a series in powers of a s (^s). Similarly, 
the more general version of the FOPT 



(33) 



is obtained by reordering (|32p in powers of a s (£so). 

Starting from (|32[) and ([3"3"]) , the new CI and FO expan- 
sions for an arbitrary scale can be obtained in a straight- 
forward way. Both the expansion coefficients and the ex- 
pansion functions appearing in the generalizations of (|20[) 
and ([50]) will depend on the scale £ . 
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5 Determination of a s (ml) 

Our objective is to compare the results of the standard 
CIPT and FOPT from [TB], [17], with the predictions of 
the new perturbation expansion presented in Section[4] To 
facilitate the comparison, we shall adopt the phenomeno- 
logical value quoted in [T7J : 

^phcn = °- 2042 ± 0.0050. (34) 

The determination of a s (m^.) then amounts to the calcu- 
lation of defined in (5) using a specific expansion of 
D{s) for s = tu^., and solving the equation 5^ — S^ en 
with respect to the coupling. 

We use the known coefficients K n from ([3|). For K 5 
we adopt the central value ([4]) with an error of ± 50%. 
Moreover, we follow the analysis made in [17] , based on 
earlier works [26], [27], [25], which leads to: 

7i = 1.21, 72 = 2.58. (35) 

The Taylor coefficients b n of the Borel transform, defined 
in (fTQ| . are 

b = 1, b t = 0.7288, b 2 = 0.6292 

6 3 = 0.7181, 64 = 0.4601. (36) 

Then the coefficients c„ appearing in the new expansion 
([T9]) , truncated at n < 4 read: 

c = 1, ci = -0.7973, c 2 = 0.4095, 

c 3 = 8.6647, c 4 = 2.2416. (37) 

The new CIPT is given by the expansion (|20p truncated 
after AT = 5 terms, with the numerical values of c„ given 
in ([37]) and the functions W„(s) defined in (j2l"]) . Inserting 
([20]) into {TJ> and using (|34|) . we obtain the prediction of 
the new CIPT: 

a s (m 2 T ) = 0.3198 ± 0.0042 oxp Ks tfrooil scale, (38) 

where the experimental error is due to the uncertainty 
quoted in (|34[) . the second error is obtained by varying 
the coefficient K$ given in ((4]) by ± 50%, and the last 
error is obtained by varying the scale £ between 1 — 0.63 
and 1 + 0.63 [TBJ. 

For the new FOPT the coefficients c„(s) defined in 
1(29]) as: 

Co (a) = 1, 

Cl ( s ) = -0.7973- 2. 6667 7 ?0 , 

c 2 (s) = 0.4095 - 2.4610 % + 3.5556 77^, 

63 (s) = 8.6647-0.6704?/o + 8.1489?7g- 3.1605 7^, 

64(a) = 2.2416- 18.1308 ?/o + 6.7302 ryg- 11.0153 r/jj, 

+ 2.107 r/o 4 , (39) 

where 770 = ln(s/so). 

The new FOPT is given by the expansion (f3T)]) trun- 
cated after 5 terms, with the numerical values of c„(s) 



given in ([39]) and VV„ defined in ([3"T]) . Inserting {30]) into 
(fTJ) and using 134|) , we obtain the prediction from the new 
FOPT: 

a s (m 2 r ) = 0.3113 ± 0.0038 cxp ± 0.0013 K , ±M scale, 

(40) 

where the errors have the same significance as in (|38[) . 

As seen from (|3"5)) and (|40[) . CIPT is more sensitive 
to the uncertainty of the last perturbative term, while 
FOPT is sensitive to the ambiguity of the renormalization 
scale. This is similar to what is obtained with the standard 
summations |16j . 

Before discussing in more detail these predictions, it is 
useful to investigate the new expansions given in Section 
|4] in the case of a physical model for the Adler function 
[T7J. 



6 Model of Beneke and Jamin 

The physical model proposed in [17] is a parametrization 
of the Borel transform B(u), consisting of one UV renor- 
malon and two IR renormalons with specified branch point 
behaviour, multiplied by polynomials. The parameters of 
the model were adjusted such as to reproduce the first 
five coefficients K n given in Eqs. ([3]) and {4]). The explicit 
expressions and the values of the parameters are given in 
Section 6 of [17], and we shall not reproduce them here. 

The new CIPT and FOPT can be constructed in a 
staightforward way: from the parameters K n of the modeQ 
one computes the Borel function {5]) truncated after a cer- 
tain number of terms N. The coefficients c„ and c„(s) at 
that order are calculated from the expansions {T9]) and 
{2"9"]) . respectively. The new expansions are given by ([2"0"]) 
and (|30p . truncated at the same number of terms A^, the 
expansion functions being defined in (|21j) and (I31j) , respec- 
tively. 

As in pPF] , we assume that the expansion of the j3 func- 
tion contains only four terms, with the coefficients given 
by ©. For a s (m1) = 0.34 the exact value of S 1 - ' , obtained 
with the PV of the Borel sum, Eq. {12]), is pi] 

6 l gl =0.2371 ±0.0060 i, (41) 

where the error is an estimate of the prescription ambigu- 
ity (cf. Eq. (6.3) of [T7J). We note that the "exact" results 
are obtained by using the contour improved Borel sum {SJ , 
with B(u) of the form specified by the model. 

The comparison of the standard CIPT and FOPT with 
the new CIPT and FOPT is seen from Figs. CD and HJ where 
we show 8^ calculated as a function of the order up to 
which the series have been summed. Fig.[T]reproduces Fig. 

4 The values from K\ to K5 are given in Eqs. ((3]l and Q; 
the next 6 values are given in Table 2 of [17] : Ke = 3275, 
K 7 = 1.88-10 4 , K s = 3.88-10 5 , K 9 = 9.19-10 5 ', K 10 = 8.37-10 7 , 
Kn = -5.19 • 10 8 , K12 = 3.38 • 10 10 ; for completeness, we 
list the next 5 parameters: K13 = —6.04 • 10 , if 14 = 2.34 ■ 
10 13 , K15 = -6.52 ■ 10 14 , K 16 = 2.42 • 10 17 , K 17 = -8.46 • 
10 17 , Kis = 3.36- 10 19 . 
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6 8 10 12 
Perturbative oder N 



Fig. 1. Values of for the model of Beneke and Jamiri 
calculated with the standard CIPT and FOPT, as a function 
of the order up to which the series have been summed. The 
horizontal band is the exact value (|41|l . 



Borel sum real part 
CIPT N=4 
CIPT N=5 
CIPT N=7 
CIPT N=10 



1.5 2 

(p (radians) 



Fig. 3. Real part of D(s) for s = m T e llfi , calculated with the 
standard CIPT, Eq. ([2]), for various numbers N of perturbative 
terms. 




6 8 10 12 
Perturbative oder N 



Fig. 2. Values of <5' ' for the model of Beneke and Jamin, 
calculated with the new CIPT and FOPT defined by Eqs. (J20J 
and (|30[) . respectively, as a function of the perturbative order 
N. The horizontal band is the exact value (|41[) . 



7 of [T7] , and shows that the standard CIPT does not ap- 
proach the true value, staying below it up to the orders 
at which the results start to exhibit large oscillations. On 
the contrary, as seen in Fig. \2k the new CIPT gives very 
good results which approach the true value with great ac- 
curacy when N increases. As concerns FOPT, the new 
approach gives results somewhat poorer than the stan- 
dard one at low orders. At large orders, when the stan- 
dard FOPT shows large oscillations, the new FOPT leads 
to values closer to the true result, but not as good as those 
obtained with the new CIPT. 

In order to understand the origin of these results, we 
calculated the Adler function D(s) for complex s along 
the integration contour. In Figs. [3]- El we present the real 
part of D(s) calculated with the standard/new CIPT and 
FOPT, for s along the upper semicircle in the definition 
Q of £(°> (s = m 2 r &V, for < <p < tt). Figs. 1 and repro- 
duce Fig. 9 of [17]. In Figs.[7]-[Tn]we present the imaginary 
part of D(s) for the same values of s. Note that the values 
along the lower semicircle follow immediately from the re- 
ality condition D{s*) = (D(s))*. By comparing Fig. [3]with 
Fig. [U one can see that, for the new CIPT, the quality of 
the approximation of the real part of D(s) improves con- 
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Fig. 4. Real part of D(s) for s = e lv , calculated with the 
new CIPT, Eq. (20j). 
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Fig. 5. Real part of D(s) for s — m^e 1 ^ calculated with the 
standard FOPT, Eq. ©. 



tinously with increasing N along the whole contour. By 
contrast, for the standard CIPT the low orders are not 
able to give a good approximation, while starting from 
N=10 the deviations increase dramatically (we can not 
show the curve for N = 15 due to these huge oscillations). 

As concerns FOPT, the comparison of Fig. [5]with Fig. 
[B] shows that the new FOPT gives a very good approx- 
imation to the real part of D(s), which improves conti- 
nously with increasing N, for if close to 7T, i.e. near the 
spacelike axis. However, the description deteriorates as <p 
approaches 0, i.e. near the timelike axis. This can be un- 
derstood by the imaginary logarithms present in the ex- 
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Fig. 6. Real part of D(s) for s = mje 1 ^ calculated with the 
new FOPT, Eq. i[30]). 
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Fig. 8. Imaginary part of D(s) for s = m^e 11 " 5 calculated with 
the new CIPT, Eq. 



pansion (|6|) of the running coupling, which are large here. 
By contrast, the approximation provided by the standard 
FOPT does not reach the same precision for ip close to it, 
i.e near the euclidian axis. 

Similar conclusions are obtained for the imaginary part 
of D(s), shown in Figs. [7]-[TDJ the new CIPT provides for 
all ip a very good approximation, which improves conti- 
nously as N increases, while the new FOPT reproduces 
well the imaginary part of D(s) in the region where the 
effect of the large imaginary logarithms is small. By con- 
trast, the standard CIPT and FOPT are not able to ap- 
proximate accurately the function at low N, and start to 
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Fig. 9. Imaginary part of D{s) for s = m^e 1 ^ calculated with 
the standard FOPT, Eq. 0. 
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Fig. 10. Imaginary part of D(s) for s = m^e llp calculated 
with the new FOPT, Eq. <[30]l. 



oscillate violently at large N. In spite of the rather poor 
local accuracy, the standard FOPT at low orders gives ac- 
ceptable results for 5^ because the region near <p = is 
suppressed by the factor w(s). 

We recall that the CI summation was proposed in [4], 
[5] in order to avoid the large imaginary logarithms, re- 
sponsible for a slow convergence of the expansion ([6]) along 
the integration countour. This slow convergence affects 
also the new FO expansion, which has as starting point 
the standard expansion ([7]). This explains the poor ap- 
proximation of the Adler function by the new FO expan- 
sion at points far from the euclidian axis. One expects of 
course an improved convergence for a smaller value of the 
coupling. This is confirmed by Fig. Ill[ where we present 
the values of 5^ calculated for a s (m^.) = 0.26 with the 
standard and new expansions. For this coupling the new 
CI and FO expansions are very close up to the 18th order, 
while the standard expansions exhibit a larger difference 
(although smaller than in the case of the physical value of 
a s (m 2 T )). 



7 Discussion 

In this Section we make several more remarks on the ex- 
pansions in powers of conformal mappings and their im- 
plications for the determination of a s . As discussed in Sec- 
tion [U the definition of the optimal variable (fT3|) requires 
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Fig. 11. Values of 5^°' for the model of Beneke and Jamin 
calculated with the CI and FO expansions, for a s (ml) = 0.26. 
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Fig. 12. Error of the determination of 8^°' as a function of 
perturbative order, for the standard and the new CI expan- 
sions. 



only the knowledge of the location of the singularities of 
B(u) in the Borel plane. The most general expansions 
based on the powers of the optimal variable are p^|) and 
(|27|) for the CI and FO summations, respectively. 

In the analysis presented in the previous Sections we 
used the modified expansions (|20|) and (|30|) . generated by 
the Borel series (fT9|) and (f29|) respectively, which explic- 
itly include the nature of the lowest (leading) renormalon 
singularities in the Borel plane. It is of interest to consider 
also the general expansions (fl6|) and (|27|) . whose respec- 
tive counterparts (fT5|) and (|26[) do not manifestly display 
this singular behaviour. According to general arguments 
|30j . we expect them to approach the true result, if the 
number of the perturbative terms is large enough. This 
is confirmed in Figs. fT2l and 1131 where we show the error 
of the determination of 5^ as a function of perturbative 
order N, for the standard CIPT and FOPT and the new 
expansions defined above. 

For small N, the expansions (fl"6|) and (f27|) give results 
similar to the standard CI and FO expansions, while the 
expansions ([20]) and (|30f approach much better the true 
result. This shows that at low N the most important ef- 
fect is the factorization of the singularities. However, for 
TV greater than 7 the effect of the conformal mapping be- 
comes visible: as illustrated in Fig. [T2l while the standard 
CI expansion starts to show large oscillations, the expan- 
sion ()16|) based on the conformal mapping of the complex 
Borel plane leads to only small oscillations around the 
true value. The improvement brought by the conformal 
mapping becomes visible rather slowly due to the strong 
leading IR singularity of the physical model considered in 

For the FO summation, the standard expansion starts 
to oscillate wildly for N > 14, while both modified expan- 
sions remain close to the true result. In this case, the ap- 
proximation provided by the expansion (|27p . which does 
not include the nature of the singularities, turns out to 
be slightly better than the modified expansion (f30|) . As 
discussed in Section [6l all the FO expansions provide a 
less satisfactory description of the Adler function near the 
timelike axis, but this region is suppressed in the integral 

©■ 

The good approximations provided by the new CIPT 
(and, at points where the convergence of the series ([6]) is 




Perturbative oder N 



Fie. 13. Error of the determination of as a function of 
perturbative order, for the standard and the new FO expan- 
sions. 



not poor, also by the new FOPT) can be explained qual- 
itatively by a theorem proved in |30j . which implies that 
the expansion (|15[) in powers of the variable w has a bet- 
ter rate of convergence than the standard expansion ([9]) at 
points of the complex Borel plane where both converge. In 
particular, the convergence is better at points on the real 
positive axis below the first IR singularity, which gives the 
dominant contribution to the Laplace-Borel integral ((8|), 
due to the exponential factor that strongly suppresses the 
contribution of large u. Therefore, the improved expansion 
of B(u) leads to a better approximation of the integral, at 
least for the perturbative orders investigated up to now. 
One might of course ask what happens if N is further in- 
creased. It can be shown that going up to the 36th order 
the new CIPT continues to be stable, while the new FOPT 
somewhat deteriorates, and this is valid also for modified 
models, having, for instance, a smaller residue of the lead- 
ing IR renormalon, compensated by an additional IR sin- 
gularity that preserve the known low-order coefficients^. 
However, a divergent behaviour at still higher orders is 
not excluded in principle, taking into account that the 
Laplace-Borel integral is performed along the cut, while 
the series (|15|) converges only at interior points. For a dis- 
cussion of this problem and explicit criteria of convergence 



We thank M. Jamin for communicating to us these results. 



I.Caprini and J. Fischer: a s from r decays 



9 



of the new expansion in the one-loop approximation of the 
coupling we refer to [19]. 

As mentioned in Section d] the inclusion of the be- 
haviour of B(u) at the first branch points is not unique. 
Of course, for a large number of terms in the expansion 
the form of these factors is irrelevant, but at low orders 
one prescription may be better than the other. In ([2"D"]) and 
(|30p . the dominant singular factors were expressed in the 
w variable. Other possible expansions will be investigated 
elsewhere. 

In order to test the choice made in this work, we use an 
idea applied in [T8] for illustrating the effect of conformal 
mappings on the convergence of power series. Consider the 
expansion (fT9|) truncated after N terms. The coefficients 
c„, n < N — 1, are calculated in a straightforward way 
from the first N coefficients K n , assumed to be known. 
When expanded as a Taylor series in the variable u, this 
expression reproduces the N coefficients K n used as in- 
put, but contains also an infinity of higher order terms, 
in particular it predicts the next coefficient K^+i- As an 
exercise, we used the expansion (|19p with 4 terms, using 
as input the coefficients K\ to K4. The coefficients cq to 
C3 are given in (|37[) . the coefficients c„ with n > 4 are 
set to zero. By reexpressing B(u) from (fT9|) as a series in 
powers of the u variable, and using ([9]) and ([TO]) , we ob- 
tain K 5 = 256, which is quite close to the value K 5 — 283 
adopted as a good estimate in [17] . So, the expansion (|19| 
is able to reproduce to a reasonable extent the higher or- 
der coefficients of the expanded function. 

This nice feature is even more striking if we go one step 
further, using as input 5 coefficients K n from ([3]) and (j4]). 
Then we have 5 nonzero coefficients c„, given in (|37p . By 
expanding B(u) from (|19|) as a series in powers of u, we 
predict the higher coefficients K 6 = 2929, K 7 = 1.73 • 10 4 , 
Kg = 3.14 • 10 5 , K 9 = 9.23 • 10 5 , K w = 6.34 • 10 7 , K n = 
-3.52 ■ 10 8 , K 12 = 2.48 • 10 10 , which agree qualitatively 
with the values of the model of Beneke and Jamin [T7] . 
given in footnote 4. These apparently miraculous predic- 
tions are explained by the fact that the new expansion in- 
corporates some features of the expanded function, which 
are known in advance. Then, with a smaller number of 
terms in the variable w, it reproduces well the higher or- 
der coefficients in the variable u. 

Finally, let us consider in more detail the predictions 
for N = 5, which is the number of terms known in the 
physical case (recall that we adopted the value of K 5 used 
in [17] )■ As seen from Fig. H for N = 5 the new CIPT 
and FOPT give results which approximate well the true 
value of 5^°' from above and below, respectively. These 
results can be understood from the plots of the real and 
imaginary part of D(s) shown in Figs. [T4TT51 The new 
CIPT and FOPT approximate very well the true func- 
tions for large ip. For intermediate values of p, CIPT and 
FOPT approach the true values comparatively well from 
opposite sides. For <p close to 0, CIPT gives definitely a 
better approximation, but this region is suppressed in the 
integral (JTJ) . This explains why the resulting values of 5^ 
given by the new CIPT and FOPT are comparable. 
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Fig. 14. Real part of D(s) for s = m^e 1¥> calculated with five 
perturbative terms in the new CIPT, Ea. pO) . and new FOPT, 
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Fig. 15. Imaginary part of D(s) for s = m^e llp calculated 
with five perturbative terms in the new CIPT, Eq. (|20p . and 
new FOPT, Eg. ([30) . 



8 Summary and conclusions 

In this paper we applied a new perturbation series for 
QCD observables, proposed in [18], to the two renormal- 
ization group summations, CI and FO, used for the ex- 
traction of a s from r decays. The new expansion, which 
replaces the standard series in powers of the coupling, ex- 
ploits in an optimal way the information about the high 
orders of perturbation theory. As discussed in detail in 
[20] , the method separates the problem of convergence of 
the perturbation series from that of its ambiguity, solved 
by choosing a prescription, which is included in the defini- 
tion the expansion functions (we adopt here the Principal 
Value). 

In the present paper we worked out in detail the CI 
and FO versions of the new perturbation expansion for 
the Adler function in massless QCD. Also, a novelty is 
the incorporation of the singular behaviour of the Borel 
transform near the first branch points by factors in the 
new variable w, as shown in ([TTJ]) and ([29]) . 

In Section [6] we illustrated the power of the new per- 
turbation theory using as an example the model of Beneke 
and Jamin [17]. As expected, the new CIPT proves to be 
superior, approaching the exact result to a very good ac- 
curacy when the perturbative order increases. The limi- 
tations of FOPT due to large imaginary logarithms along 
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the integration contour in the complex plane are clearly 
illustrated. 

From the predictions ([38]) and (j40|) . adding an uncer- 
tainty of of 0.003 due to the power corrections [17] , we 
obtain: 

a s (m 2 T ) = 0.3198 ±°;°Jg > new CIPT, 
a.(m£) = 0.3113 ±°;°^, newFOPT. (42) 

As discussed in Section [5] the dominant contribution to 
the error is due to the uncertainty of the last perturbative 
term in the case of CIPT, and to the ambiguity of the 
renormalization scale in the case of FOPT. 

It is remarkable that the difference between the cen- 
tral values in (|42|) is only 0.009, while for the standard 
expansion the difference is 0.024. The new expansions re- 
move thus the most intriguing theoretical discrepancy in 
the determination of a s from r decays. We note that both 
values in ([4"2"]) are closer to the standard FOPT than to 
the standard CIPT. So, our analysis indirectly confirms 
the criticism of the standard CIPT made in [17]: although 
the renormalisation group summation is more accurate, 
the low order perturbative terms are not able to describe 
the high order features of the series. The new CIPT brings 
an improvement precisely at this point. 

According to the last remark made in Section[7] for the 
model considered in [17] the new CIPT and FOPT with 
N = 5 terms give comparable predictions, which approx- 
imate the true result from opposite sides. If this model 
describes correctly the physical situation, then the true 
result is expected to be between the two predictions in 
([42]) . With this assumption, one may take the weighted 
average of these two values, which leads to 

a s (m 2 T ) = 0.315 ±g;ggl , average. (43) 

However, in order to avoid any bias related to a specific 
model, we take as best result the value given in ([42]) by 
the new CIPT: 

a s (m 2 T ) = 0.320 l&SSJ. (44) 

We recall that this prediction is based on the new contour- 
improved expansion defined in Eq. (|20p . This expansion 
reproduces the known perturbative coefficients K n order 
by order, includes the information about the dominant 
singularities of the Adler function in the Borel plane, and 
is based on an optimal expansion of the Borel transform, 
which converges in the whole complex it-plane up to the 
cuts along the real axis. 

The result ([44]) coincides practically with that obtained 
in [17] using the standard FOPT. However, our prediction 
has a more solid theoretical basis, being free of fortuitous 
compensations of terms related to large imaginary loga- 
rithms, like in FOPT. Also, it is based on a systematic 
perturbation theory, and its uncertainty is related mainly 
to the error of the last perturbative term. So, the accu- 
racy of the prediction is expected to increase when more 
perturbative terms for the Adler function in QCD will be 
available. 
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